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Abstract 



> 

!\^ ■ We study a generalized ladder resummation in the superfluid phase of 

^^ . the nuclear matter. The approach is based on a conserving generalization of 

C**** I the usual T-matrix approximation including also anomalous self-energies and 

propagators. The approximation here discussed is a generalization of the usual 

mean-field BCS approach and of the in medium T-matrix approximation in 
en \ the normal phase. The numerical results in this work are obtained in the quasi- 

particle approximation. Properties of the resulting self-energy, superfluid gap 

and spectral functions are studied. 
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1 Introduction 



One of the most general properties of fermion systems with attractive interactions 
is the transition to a superfluid state at finite density and low temperature. It is 
expected that such a phenomenom happens also for the strongly interacting system 
of nucleons, i.e for cold nuclear matter. Calculations based on nucleon-nucleon 
interaction predict very large values of the superfluid gap. Typically around 5 — 
lOMeV for the isospin singlet (T = 0) partial- wave 3 Si - 3 D 1 0, g g, |. The actual 
value of the superfluid gap is a matter of debate, because of expected screening 
and self-energy corrections. One has to note that also some of the phenomenological 
potentials fitted to the pairing properties of finite nuclei give significant values of the 
superfluid gap in nuclear matter 0. The study of T = pairing in the systematics 
of TV ~ Z nuclei became possible with advent of radioactive beam facilities. This 
lead to a resurgence of the study of the nuclear mass systematics || |7], ||. Thus 
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it is of importance to obtain results on the nature of the nuclear pairing (T = 
or T = 1) and the value of the superfluid gap for symmetric nuclear matter at 
saturation density and below. It appears that the mean-field gap equation without 



medium modifications is unrealistic || [TO], [TT], [12], [13], |14j and hence the best strategy 
would be to calculate a density dependent gap for nuclear matter (including relevant 
many-body corrections) and use it in a local density approximation for calculations 
in finite systems H |7|. 

The neutron rich nuclear matter is used in modeling of the crust and of the core 
of neutron stars. It is generally believed that such an asymmetric nuclear matter 
is superfluid, with different kinds of superfluid gap appearing in the vast range of 
densities present in the neutron star. The value of the superfluid gap is of importance 
for fundamental problems in the neutron stars, the formation of glitches, the value 



of the viscosity and the cooling rates in different scenarios |15| . Again one of the 
possible approaches is to calculate the superfluid gap from bare NN interaction 
using the Brueckner-Hartree-Fock approximation to get single-particle energies for 



the gap equation |T6], [17], [TT], [I3 [L9[ [13|, [14], |2(J . However, in order to obtain reliable 



estimates for the superfluid gap in the neutron matter we need to have under control 
in medium many-body corrections to the gap equation. 

As noted above strong modifications of the mean-field result for the superfluid 
gap are expected in the nuclear matter due to screening effects || pi], |ll|], modi- 
fications of the effective mass and self-energy corrections [|T^, [TJ], [Tj|. One of the 



motivations of the present study is to investigate another source of in medium cor- 
rection to the value of the gap. These corrections occur in a generalization of the 
resummation of ladder diagrams to include also anomalous propagators. In the fol- 
lowing we study the formalism and give numerical estimates of the corrections to 
the anomalous self-energy These additional terms introduce an energy dependence 
in the superfluid gap, but the modifications of the value of the gap at the Fermi 
surface are not as dramatic as from the other medium effects. 

The correction to the binding energy due to the superfluid rearrangement of the 



ground state is believed to be small |21|] . However, some of the calculations using 
realistic nuclear forces predict quite large values of the superfluid gap in the nuclear 
matter. A large superfluid gap could lead to modifications of the normal part of 
the self-energy and the spectral function. In the following we study a consistent 
approximation treating on equal footing the normal and anomalous part of the self- 
energy We find that important modifications of the single-particle and two-body 
propagators appear if the superfluid gap is large. These significant modifications in 
the superfluid present, a posteriori, an important reason to consider the generalized 
formalism discussed in this work. Also the expansion of the ground state energy or 
other quantities around the wrong ground state is not satisfactory for a theory aiming 
at the description of the many-body problem from first-principles, using free nucleon- 
nucleon potentials. The incorrect ground state could also lead to instabilities in 



the actual iterative numerical solution of the many-body equations, related to the 
appearance of the Cooper instability |22|] . 

Mean-field approaches give a qualitatively correct description of the formation 
of the superfluid gap by the BCS mechanism, but fail in the resummation of the 
hard core. Recently superfluid nuclear matter was studied in an approach starting 
from the in medium T-matrix approximation P%| . The resummation of ladder di- 
agrams in the T-matrix can be used to deal with the hard core in the interaction 
potential. The T-matrix approximation for the self-energy was studied intensively 
in the last decade in normal nuclear matter. It is different from the usual G-matrix 
approximation. The T-matrix formalism, also called self-consistent Green's func- 
tion approach, can be used directly at high temperature, above the superfluid phase 
transition [F2J], |2j| ^B[ [2TJ. It allows to study in a self-consistent way the one-particle 
self-energies and spectral functions [|4], |2j| [2j| ^9| |30| , the two-particle properties 



and in medium cross-sections p3L ER 12b], E7L KJ], |2j|, the onset of superfluidity 



32], ^, |27], |23], |^] , the self-energy corrections to the superfluid gap |12| . The treat- 
ment of the old question of saturation properties of nuclear matter in the T-matrix 
approximation is at the present stage not superior to the most recent G matrix or 
variational calculations, including realistic interactions, three-body and three-body 
forces corrections. However, the T-matrix self-energy leads to reliable results for 
the single particle properties, in particular it gives a consistent value for the Fermi 
energy fulfilling the Hugenholz-Van Hove theorem [3~3, P7| . 



Besides the above motivations to develop the T-matrix approach for nuclear 
matter, this approximation seems to be the most natural starting point for the 
study of the superfluid phase of the nuclear matter. The appearance of a singularity 
in the T-matrix at zero total momentum of the pair and at twice the Fermi energy 
signals the formation of a long range order and defines the Thouless criterion for the 
critical temperature. It is in fact equivalent to the critical temperature corresponding 
to the appearance of a nontrivial solution of the BCS gap equation. The same is 
true if off-shell propagators are used in the T-matrix ladders. It corresponds to a 
generalization of the gap equation to one using full spectral functions, including the 



imaginary part of the normal self-energy |12j. The study of the singularity in the 



T-matrix equations at finite temperature, in the normal phase, allows to identify 
the critical temperature and the precritical modifications of the spectral function, 
in medium cross-sections and density of states [321 [33], ^, |35|, [2D], |2j| |23L [36] . 



The generalization of the T-matrix approximation to the superfluid phase ap- 
pearing below T c was discussed in Ref |[23|| . The approach is based on the observation 



that the formation of the superfluid order parameter requires a long range order. 
This long range order, representing the formation of Cooper pairs at zero momen- 
tum and twice the Fermi energy, corresponds to a singularity in the T-matrix for 
the same energy and momentum as in the Thouless criterion. Thus also below T c we 
expect the singularity in the two-body propagator at twice the Fermi energy. The 



kernel of the T-matrix equation is modified in the superffuid so that the singularity 
of the T-matrix equation is again equivalent to the BCS gap equation for non-zero 
values of the order parameter [f23] . 



In the present work we investigate a different approach to the T-matrix resum- 
mation in the superffuid. It is a generalization of the ordinary T-matrix ladder 
diagrams which includes also anomalous propagators. Thus normal and anomalous 
self-energies are calculated in a unified way. The approximation deals at the same 
time with ladder diagrams resummation for the self-energy and with the appearance 
of the order parameter. The properties of the generalized T-matrix, self-energy and 
spectral function are discussed. A justification of the heuristic procedure of Ref. 
23] is then obtained if the superffuid order parameter is restricted to the BCS con- 
tribution. The practical calculations are more difficult in the generalized scheme 
here discussed. The number of propagators and T-matrix components is doubled 
because of the appearance of the off-diagonal, anomalous propagators. The numer- 
ical results presented in this work are obtained in the quasi-particle approximation, 
starting from mean-field BCS propagators. 

2 Green's functions in the superfluid 

2.1 Notation and formulas for the normal phase 

We consider infinite homogenous nuclear matter interacting through a two-body 
potential. The energies are defined with respect to the chemical potential \i 

a',/3',a,l3 

In the real time formalism the Green's functions are defined on a contour in the time 



plane 37 . It means that as functions of a single valued time variable the Green's 



functions acquire a matrix structure 

Q ( x t'X t ) = ( ^a/^ 1 '* 1 '^ 2 '* 2 ) G a/3 (x 1 ,t 1 ;X 2 ,t 2 )\ ,£) 

\G a p(Xi,ti',X2,t2) Gap[Xl, t\\ X2, £2) / 

Where G^g(xi,ti]X2,t2) = i(T'9 a (xi,ti)Wa(x2,t2)} is the chronological Green's 
function, G a is the anti-chronological Green's function and GS(xi,£i;x2,£2) = 

i($ a (xi, ti)^fp(x2, t 2 )) and G< 3 (xi,tx; x 2 , t 2 ) = -i(^ g (x 2 , t 2 )^ a (xi, ti)) are the cor- 
relation functions. In a homogeneous system, after Fourier transforming in the rela- 
tive coordinate and time, single-particle Green's functions depend on a single energy 



and momentum. Noting that the Green's function are diagonal in the spin-isospin 
indices, we can write 

G a/3 (p,u}) = 5 af3 G(p,uj) . (3) 

A single scalar function (in spin-isospin indices) G is sufficient to describe a spin 
and isospin symmetric nuclear matter. The scalar correlation functions are written 
using the spectral function 

G<(p,u)=if(u)A(j>,u) (4) 

G>(p,u>) = -i(l-f{u))A(p,u) . (5) 

f(u)) = 1+ex lg( w _ ) is the Fermi distribution at temperature T = 1//5. The spectral 
function A can be related to the retarded Green's function 

A{p,u) = -2ImG + (p,uj) 

= i{G > {p,u)-G<{p,u)) (6) 

-2Im£+(p,u;) 



(u - p 2 /2m - ReS+(p, u) + fi) 2 + Im£+(p, uj) 2 



(7) 



where £ + (p, u>) is the retarded self-energy. 

The T-matrix equation in the normal phase is 



\Pl T 5/3' a /3( P '^)IP) = K*'/3'a/3(P> P 



+ ^Iwj-Jwf v ^'^ k) 



d 3 k r d 3 q 

I ; 

(k|^ ± (P,a;)|q)(q|T± a/3 (P,a;)|p') , (8) 

where Q is the disconnected retarded two-particle propagator 

du r du 



(p|^(P^)|p) = (2 7 rm P -p)/^/ 



2vr 
^G<(P/2 + p, u" - u/)G<(P/2 - p, J) 

-G > (P/2 + p, J' - lj')G > {P/2 - p, J))/(uj - J' ± ze) 

= (27r) 3 5 3 (p-p')^ ± (P,u;,p). (9) 

Note that it is sufficient to solve a single equation for the retarded T-matrix instead 
of a matrix equation in the indices on the time contour in the complex time plane 

©■ 

A full structure in the spin-isospin indices and relative angles of momenta must 

be kept. However, usually a partial wave expansion of the T-matrix equation is 



performed 

i, J Z7T 

g ± (P,u,k)(k\T^ l ST)± (P,u;)\p) (10) 

after angle averaging the Kernel Q 

g ± (p,co,p) = J^g ± (p,co, P ). (ii) 



2.2 Anomalous Green's function 

In the superffuid phase the ground state of the system has a nonzero order parameter, 
corresponding to bound Cooper pairs PS[. This leads to new Green's functions 



p ( T +. r ^.(Kp^iM^^M) F< /3 (x 1 ,t 1 ;x 2 ,t 2 )\ n9 x 

* a 0{.X 1 ,t 1 ,X 2 ,t2) - I p > / +.„ +\ paf T y. . _ y. \ J • l^J 

V ^a/SV^l) H) ^2, £2 J -r a/3^1, tlj #2, t 2 ) / 

F c is the time ordered anomalous Green's function 

F^(xi, h; , x 2 , h) = HT^JxxM^fa* h))) 
and F a the anti-chronological one and where 

F< p (x u ti;,x 2 ,t 2 ) = i(y a (xuti)-$0(x 2 ,t 2 )) (13) 

F> p (x u ti;,x 2: t 2 ) = -i(y f3 (x2,t 2 )y a (x 1 ,t 1 )) . (14) 

We consider a homogeneous infinite system so that the anomalous Green's functions 
depend on a single momentum and energy. Analogously as for the normal Green's 
functions we can write 

F^fau) = if(u)B af) fau) (15) 

F> p (p,u) = -i{l-f{u))B a p(p,u) . (16) 

The anomalous spectral function is related to the imaginary part of the retarded 
propagator 

B aP (jp,u) = -Im(F+ f} (p,oj + ie)-F+p(p,u;-ie)) (17) 

= i{F%,(p,u>) - F$(p,u)) (18) 

It should be noted that the spectral function for the diagonal part of the Green's 
function is modified in the presence of the off-diagonal self-energy A. We shall 
denote it by 

A s (p, u) = -21mG + {p, u) (19) 
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and reserve the notation A(p, u>) to the spectral function obtained by putting A = 
(Eq. 0). 

The spin-isospin structure of the anomalous Green's function is assumed to be 
of the spin (isospin) singlet or triplet kind. We write 

F af3 (p,Lj)=A af3 F(p,u;) . (20) 

In general the matrix A a/ g could depend on the momentum, energy and relative 
directions of spin, isospin and momentum. To simplify we use in the following angle- 
averaged double propagators in the T-matrix ladder. The matrix A in spin-isospin 
fulfills 

A^A^ = \D\ 2 5 a , (21) 

for time-reversal invariant states. Without loss of generality we can put \D\ 2 = 1. 
More specifically the spin (isospin) part of A a/ 3 is 

d x -\- lUy icl z \ (00\ 

for the triplet gap (d 2 , + d 2 + d 2 z = 1) and 

-, o) < 23 > 

for the singlet one. Together with the choice of a diagonal normal Green's function 
we describe the propagators using scalar functions in spin-isospin indices F <> (p,u) 
and G <> (p,u) or the corresponding retarded propagators and spectral functions. 

2.3 Ladder resummation in the superfluid 

Ladder resummations in the superfluid have been considered in the description of 
high T c superconductors ||39|| . A thermodynamically consistent scheme which reduces 



to the usuall T-matrix equation above T c can be constructed. The simplest way is 
to introduce a generalised T-matrix with aditional indices indicating if the incoming 
line is anomaleous or normal, following the Nambu formalism for superconductors. 
Reset icting oneself to two normal or two anomalous propagators in the ladder leads 
to the following additional matrix structure in the T-matrix 

l<p|^(P,-n)|p> <p|T^(p,-fi)|p')J 

(V{p,p') a >p> a p \ t d 3 q v^fyfoqWyf 

I V(p, v') al0a , ) + J (2tt)3 ^ I V(p, q) a ^ W5 ' 

8 w 6 ss ,g±(V,n,q) W± 7 , y (P,Q,q) 
±t # / 

7^7 8 



n%, 5 , (p, -n, q ) 5^5 S5l g±(p, -n, q ) 



where 



<q|r^(P,n)|p> ( q \L± 5 ,jPMp) \ 

( q \L±4 af) (P,-n)\p) ( q \T^,jP,-Q)\p')) [ ] 



du f duo' 



*SJ?,n,P) = f%f 



2n 



-F>,(P/2 



CTf\ 



(F< (P/2 + p, J' - oj')F< (P/2 - P , u/) 
p,u/' - u}')F£ s (P/2 - p,u'))/(tt - uj" ± ze 



(25) 



denotes two disconnected anomalous propagators. L^ is the off-diagonal part of the 
generalized retarded T-matrix (Fig. |l|). 

The anomalous part of the ladder Ti a p^$ does not mix different values of the total 
spin (isospin) of the pair. It can be seen by writing the matrix structure of H a p in 
the basis of the total spin (isospin) of the pair and of its third component. 



/ 



n 



(d x + d y ) 2 -i\/2d z (d x + idy) 
i\^2d x (d x + idy) d 2 x + d 2 y — d\ 
—d? z —id z [d x — idy) 



-dl \ 

—i\[2d z [d x — idy) 
(d x - idy) 2 
D 2 ) 


H (26) 


et pairing with the matrix in the cor 


nponents 1 + , 1°, 1 _ , 0° of the total 



spin (isospin) of the pair and 



n 



/ 


-1 


°\ 


1 








-1 








V o o 





1/ 



H 



(27) 



for the singlet pairing. Using the structure of the anomalous Green's function (Eq. 
20[ ) and assuming a dependence of the angle averaged spin-isospin structure (the 
matrices in Eqs. 26, 27) only on the value of total momentum one finds for central 
forces 

£ H a ,p. ld x n i5aP = Sjfarf x H , (28) 

where x stands for momentum integrals with the interaction potential or the T— matrix 
as occurring in the iteration of the Eq. (0). Thus partial wave decomposition can be 
applied after angular averaging of the intermediate uncorrelated two normal or two 
anomalous propagators in the ladder. The resulting generalized T-matrix equation 
has a matrix structure corresponding to (Eq. |24|) in each partial wave 

/ ( P \T± ST) (p,n)\p) (p\Lf JST) (p,n)\ P ) \ = 

{ (p\Lf JST) (P, Q)\p) (p\T±; sT) (P, -Q)\p) ) ~ 



V(j S T)(p,p) \ | fq 2 dq /V {J sT){p,q) 

V(jsT)(p,p')J J 2vr 2 V V( JST )(p,q) 

( g ± (p,n, q ) 7^(p,n,g) 

/(g|T± 5T) (P,fi)b') (q\Lf JST) (P,Q)\p) \ 
{ (q\Lf JST) (P,Q)\p) (q\T$ ST) (P, -Q)\p) ) 



(29) 



We used also (p|L ± (P, £l)\q) = (p|-^ ± ^(P,—£l)\q). In the following we skip the 
spin-isospin indices in the equations assuming a partial wave decomposition in the 
generalized T-matrix. 

2.4 Self-energy 

The self-energy can be defined generalizing the self-energy in the T-matrix approxi- 
mation for the normal phase. The diagonal part of the self-energy is defined by the 
T-matrix approximation to the two-particle Green's function (Fig. |2|) 

S<>(p^)=./^^((p-k)/2|T<>(^ + ^,P + k)|(p-k)/2)G<>(A;^') . 

(30) 
Analogously we can define the T-matrix part of the off-diagonal self-energy (Fig. §) 

A<>(p^)=./^^((p-k)/2|L<>(c, + c,',p + k)|(p-k)/2)P<>(fc^') . 

(31) 
The above definition of self-energy is a ^-derivable approximation |4(| (Fig. |]). It 
leads to thermodynamically consistent results. In Eqs. (|3T]) and (p0|) the subscript a 
denotes the antisymmetrization and we have not written explicitely the spin-isospin 



indices. It can be checked that the formula ([31]) conserves the singlet (triplet) 
structure of the order parameter given by F. The approximation ([31]) does not 
include the usual BCS part of the anomalous self-energy 

^bcs(p) = ~*J^J (03^' k ) F< ^ <") " ( 32 ) 

This term can be added explicitely to the the off-diagonal part of the self-energy 
(Fig. H). Moreover the $ derivability of the approximation is conserved. As it will 
turn out the BCS part of the anomalous self-energy is dominant, and in the first 
approximation one can neglect the two-body contribution (|31~D to the superfluid 
gap. However the approximation scheme consisting in keeping the T-matrix form 
of the diagonal self-energy and the BCS form of the off-diagonal self-energy is not 
^-derivable. 



2.5 Quasi-particle approximation 



The full solution of the set of equations (|24|) for the generalized T-matrix, the 
normal self-energy and the superffuid gap requires a self-consistent iterative solution. 
Below we present results of a simpler calculation. It starts with the mean field 
approximation for the normal self-energy (Hartree-Fock approximation) and the 
BCS approximation for the superffuid gap 

^bcs(p) =-J j^y 3 V(p, k) — A BCS (k) . (33) 

The starting spectral function for the calculation of the T-matrix is 

A BCS (p, co) = 2n (%^*{u - E p ) + ^g^% + T p )) (34) 



where £ p = p 2 /2m+T lHF (p)—fi and E p = J& + A p . The corresponding off-diagonal 
spectral function is 

B(p, cu) = -2vr^ (5(u - E p ) - 6(uj + E p )) . (35) 

Zh, p 

The generalized T-matrix equation is solved for the S"-wave Yamaguchi interac- 
tion ETJ. This oversimplified interaction is used to reduce the numerical difficulties 



involved in the solution of T-matrix equations and in the calculation of the self- 
energies. First the imaginary parts of the diagonal and off-diagonal self-energies are 
calculated 

l^ M = J^^(f(E p) + ^ + E p )) 

((p - k)/2|ImT(p + k,cu + E p )\(p - k)/2) A 
+^^ E (f(-E p ) + b{u - E p )) <(p - k)/2|ImT( P + k,u - E p )\(p - k)/2} A ) (36) 

and similarly for IirAtm (b(u) is the Bose distribution). Then from the disper- 
sion relations the real part of the dispersive contribution to the self-energy can be 
obtained 

/ Re£(p, u) disp \ = f du / ImS(p, u') disp \ 1 , . 

\ReA(p,uj) TM ) ~ J 2tt \ImA(p,w) TM ) u-u' ' l ' 

To the dispersive part of the self-energies one has to add the mean-field self-energy 
Ti HF and the mean field superffuid gap A BC s 

ReS(p, u) = ReS(p, u) disp + Y>hf{p) 
ReA(p, u) = ReA(p, u) TM + A BCS (p) (38) 

that have been calculated already to obtain the spectral functions (|34]), (|35|). 
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3 Results for the T-matrix 

3.1 Energy gap in the T-matrix and the two-particle spec- 
tral function 

The equations for the generalized T-matrix are solved with the quasi-particle ansatz 
for the spectral functions A and B. We present results at the temperature of 3MeV 
with the BCS gap of lO.lMeV. The results do not change appreciably when reducing 
the temperature further. In Fig |5] is shown the imaginary part of the diagonal and 
off-diagonal parts of the generalized T-matrix. From the BCS ansatz for the spectral 
functions follows a gap in two-particle excitations around the the Fermi energy. For 
zero total momentum of the pair, this forbidden region is twice the superfluid gap. 
At the edge of the two-particle gap the generalized T-matrix has a singularity as 
function of energy, similar to the singularities in the one-particle spectral function 
in the BCS approximation. 

In fact the imaginary part of the T-matrix is related to the two-particle propa- 
gator in the T-matrix approximation. The one-time two-particle spectral function 
is 

(p\A 2 (P,Q)\p) (p\B 2 (P,n)\p) 

( P \B 2 (p,-n)\p) ( P \A 2 (p,-n)\p) 

J (27r) 3 i (2vr) 3 V V~^[p,p , 

(q\T±(P,n)\k) (q\L±(P,n)\k) \ (V~\k,p) 

(g|T ± (P,-fi)|A;} {q\T±\P, -Q)\k) J \ V^fcp 



(39) 



The spin-isospin indices are omitted, and the two-particle spectral function can be 

projected on a definite total spin or isospin of the pair. Also the relative momentum 

can be projected on states with definite angular momentum. In Fig. |B| we present the 

two particle spectral function projected on partial waves occurring in our separable 

interaction 

d 3 kd 3 p 



(M P- H)) = / —rfg(k)g(p)(k\A 2 (P, Q)\p) , (40) 



where g(k) = l/(k +/3 ) is the formfactor of the Yamaguchi interaction [41]. We use 
the same projection for the off-diagonal two-particle spectral function B 2 in different 
channels and for the uncorrelated BCS two-particle propagators. 

The energy gap present in the two-particle BCS propagator is visible also in the 
T-matrix approximation (Fig. |^). It corresponds to a minimal energy to excite a 
two-particle pair of twice the single-particle energy gap. A similar gap appears in 
the two-particle anomalous propagator (Fig. |^) and for nonzero momentum of the 
pair (Figs. |, |). 
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3.2 Singularity in the T-matrix 

The T-matrix is singular for the total energy of the pair equal to twice the Fermi 
energy and zero total momentum. It is a generalization of the Thouless criterion 
for the critical temperature to the superfluid phase. Indeed the imaginary part of 
T and L is always zero at Q = and P = and the inverse of the real part of the 
generalized T-matrix 

Mk'\( T r £) W^ ( 41 ) 

\-Ll 1 J p =on=Q 

has a zero eigenvalue 

<^''(I rU^M-AmH' (42) 

A(fc) (with spin-isospin indices omitted) is a solution of the mean-field gap equation 
. , . f d 3 k Ir . ,.du)du A(u, k)A a (u, k) /_ „, x ,, a>,„, „ ,.„, 

A(rt - / (2^< * W w w t 1 - /(w) - /(w ) )AW = ° (43) 

if A s is the full spectral function with only the BCS contribution to the off-diagonal 
self-energy Abcs- in particular the inverse generalized T-matrix with mean- field 
quasi-particle propagators (|51D, d55|) has a zero eigenvalue corresponding to the 
BCS mean-field superfluid gap (|33|). In Fig. [K] is shown the real part of the 
inverse determinant of the generalized T-matrix for the two S partial waves. Clearly 
the superfluid gap is formed in the 3 Si channel (for the chosen interaction). The 
generalized T-matrix in this channel shows a singularity at zero total momentum 
and at twice the Fermi energy for all temperatures below T c . 

4 Self-energies and spectral functions 

4.1 Single-particle energies 

The position of the single particle pole uo p when including only the real part of the 
self-energy can be obtained from the solution of 

p 2 

u p = h ReS(p, Up) - fj, . (44) 

2m 

In Fig. |TT] it is compared to the position of the quasi-particle pole £ p when including 
only the Hartree-Fock self-energy There is a significant difference between the mean- 
field single-particle energy and the single-particle energy (44) for low-momenta, due 



to the dispersive part of the self-energy as obtained in the T-matrix approximation. 
It is a reflection of the presence of short range correlations for nuclear interactions. 
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This modification of single-particle energies below the Fermi energy leads to impor- 
tant corrections to the binding energy per particle in nuclear matter. A resummation 
of ladder diagrams in the form of the G-matrix (or the T-matrix) is necessary to 
obtain reliable results for the binding energy |42|, |21| . 

In the superfluid the presence of the off-diagonal self-energy A(p, uj) leads to the 
splitting of quasi-particle poles in the spectral function A s (p,u) 



E p = ±^/(£ p + Re£(p, E p )f + | A(p, E p ) | 2 . (45) 

As shown in Fig. |TI| the superfluid quasi-particles, present a gap in excitations 
around the Fermi momentum. Far from from the Fermi momentum the domi- 
nant quasi-particle pole of A s (p, uj) approaches the quasi-particle pole uj p of A(p, uj). 
In Fig. [12] the position of the quasi-particle poles obtained including the full off- 
diagonal self-energy (Eq. ^5|) is compared to one where only the BCS part of the 
superfluid gap is taken 



E p = ±yl(£ p + ReS(p, E P )Y + | A BCS (p)\ 2 • (46) 

The two energies are very close. At the scale of Fig. |11] they cannot be distinguished 
and only close to the Fermi momentum a small difference can be seen in Fig. [12|. It 
is due to the fact that the dispersive part of the off-diagonal self-energy A TM (p, uj) 
introduces only a small correction to the mean-field superfluid gap A BC s(p)- 

4.2 Imaginary part of the self-energy 

In Fig. ^ we plot the energy dependence of — ImE + (p, uj) for several values of the 
momentum p. A characteristic feature is the very strong reduction of the damping 
in the region around the Fermi energy, already at T = 3MeV. It is related to the 
appearance of an energy gap for excitations of twice the value of the superfluid gap. 
On top of that the usual reduction of damping due to the restriction of the phase 
space appears, leading to the formation of quasi-particles around the Fermi surface. 
The quasi-particle nature of the excitations can be judged by plotting the single 
particle width 

T(uj, P ) = -2lmS(p,uj) (47) 

at the quasi-particle pole. The single particle width is small around the Fermi energy 
for the quasi-particle pole uj = oj p and the poles u = ±E p (Fig. [b|). The formation 
of the energy gap for excitation of particle pairs and the reduction of the phase 
space around the Fermi surface leads to the appearance of sharp quasi-particles for 
momenta close to the Fermi momentum. 
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4.3 Superfluid gap 

The superfluid gap acquires a contribution Atm from the T-matrix diagram. The 



real part of the T-matrix gap is obtained from the dispersion relation (^3). In Fig. [15 
is shown ReA^Af (p, k>) for several values of momenta. Its value is small for energies 
corresponding to the Fermi energy and for momenta close to the Fermi momentum. 

In Fig. [TB| is shown the imaginary part of the T-matrix contribution to the 
superfluid gap. The imaginary part of the dispersive contribution to the superfluid 
gap IthAtm{p,^) shows a gap around the Fermi energy similarly as the imaginary 
part of the diagonal self-energy This means that the imaginary part is not modifying 
the value of the superfluid gap around the Fermi surface. 

In Fig. [I7| the mean-field superfluid gap A BC sij>) is compared to the superfluid 
gap containing also the T-matrix contribution at the quasi-particle poles uj p or E p . 
Close to the Fermi momentum the mean-field value of the superfluid gap is not mod- 
ified significantly. Only for small momenta the T-matrix contribution to A(p,E p ) 
shows up. 

4.4 Spectral functions 

The quasi-particle nature of the excitations around the Fermi momentum can be 
judged from the spectral functions A(p,u>) and A a (p,u). In the upper panel of Fig. 
|18| are shown the spectral functions for p = (below the Fermi surface). The main 
strength of the spectral function is concentrated around the quasi-particle peak at 
uj = — 56MeV = uj p ~ — E p . However for this momentum the spectral function 
is relatively spread. The quasi-particle approximation would be of limited validity. 
Above the Fermi energy the strength of the spectral function is small. It is slightly 
larger for A s than for A because of the contribution of the second quasi-particle pole 
of the full spectral function A s (p,uj) in the superfluid. 

The situation is reversed in the lower panel of Fig. |TH| where the spectral functions 
are plotted for the momentum p = 350MeV (above the Fermi momentum). The 
quasi-particle pole of A and the dominant quasi-particle pole of A s are located above 
the Fermi energy. Below the Fermi energy there is a small contribution from the 
background strength of the spectral function to A and A s and a small contribution 
from the second pole to the full spectral function A s . 

In the middle panel of Fig [18| are shown the spectral functions for a momentum 
close to the Fermi momentum. Very sharp quasi-particle peaks are visible for A(p, uj) 
at uj = u p and for A s (p, uj) at uj = ±E p . For this momentum the difference between 
A and A s is the most pronounced, because uj p is different from E p and both poles 
of A s have comparable strength. 

In Fig. [19] is plotted the off-diagonal spectral function B(p,uj). The anomalous 
spectral function is an odd function of energy. B(p, uj) has two quasi-particle poles on 
both sides of the Fermi energy. The quasi-particle poles are very sharp for momenta 
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close to the Fermi momentum. Only for small momenta the off-diagonal spectral 
function shows two relatively broad peaks. 

5 Conclusions 

We present an approach which allows for the resummation of ladder diagrams in 
the form of the in-medium T-matrix and the treatment of the superfluid phase of 
the nuclear matter at the same time. The approximation presented in Sect. |2.3| is a 
generalization of the usual T-matrix resummation in medium to temperatures be- 
low T c . In the superfluid phase the Green's functions, self-energies and the T-matrix 
acquire additional indices corresponding to anomalous propagators. The approach 
is ^-derivable and hence it is thermodynamically consistent, like the self-consistent 
T-matrix approximation. The partial wave expansion can be performed for central 
forces. The corresponding self-energy is calculated assuming the T-matrix approxi- 
mation for the two-particle propagator and the two-particle anomalous propagator. 
The off-diagonal self-energy can be supplemented with the mean-field BCS con- 
tribution which turns out to be dominant. The addition of the BCS gap to the 
off-diagonal self-energy does not spoil the $ derivability of the approximation. The 
spin-isospin structure of the T-matrix part of the superfluid order parameter is the 
same as its mean- field part. 

In this exploratory work we use a simple separable interaction to illustrate the 
method by numerical results. The calculations are performed in the approximation 
where the normal and anomalous propagators in the generalized T-matrix ladder 
and in the self-energy diagrams are of the BCS mean-field form. It represents the 
first iteration in the calculation of the self-consistent set of equations with off-shell 
normal and anomalous propagators. It turns out that the real part of the diagonal 
self-energy changes in the first iteration significantly from its mean-field form. On 
the other hand the value of the superfluid gap is only slightly modified by the 
addition of the T-matrix contribution in the first iteration. The imaginary part of the 
diagonal self-energy is reduced around the Fermi energy due to the appearance of an 
energy gap for excitation of twice the value of the superfluid gap. The quasi-particle 
approximation is a reasonable description of excitations close to the Fermi surface, 
with single-particle energies modified by the resummation of ladder diagrams. 

Let us comment on the influence of the use of more realistic forces on the results. 
The importance of the ladder diagram resummation for the real-part of the self- 



energy is of course important for any realistic interactions [ i2| , pTT , |35, |36], |26 |. The 
value of the superfluid gap is expected to be much smaller for realistic nuclear 
potentials at normal nuclear density. The validity of the partial wave expansion for 
non-central forces in the superfluid remains to be checked. However since we expect 
that the dominant part of the off-diagonal self-energy comes from BCS diagram the 
partial wave expansion in the generalized T-matrix calculation would not introduce 
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important modifications for final self-energies. 

Another open question is the effect of the self-consistency of the calculation on 
the spectral functions and self-energies. The equations for the self-energies and 
the generalized T-matrix should be iterated instead of using only the quasi-particle 
mean-field propagators. In the self-consistent iteration the superfluid gap is reduced 



12] and the energy gap in the superfluid is no longer sharp. The imaginary part of 



the self-energy will also be reduced around the Fermi energy but not as sharply as in 



Fig. |14[ We note that only the use of self-consistent propagators and self-energies in 
the generalized T-matrix approximation guarantees the thermodynamic consistency 
of the results ||27||. As a result of the self-consistency in the off-diagonal self-energy 



the generalized T-matrix will not have a singularity at twice the Fermi energy and 
zero total momentum. As described in Sect. |3.2| this singularity is present only if 
mean-field BCS off-diagonal self-energy is used in propagators of the ladder. The 
inclusion of two-particle contribution to the order parameter modifies this criterion 
for long range order in the two-particle propagator. 



Finally comparing to results of Ref. [23| we found a similar behavior of the real 



part of the generalized T-matrix if using BCS propagators, i.e. a singularity at twice 
the Fermi energy and zero total momentum. However the imaginary part of the T- 



matrix is very different form the results in Ref. |[23|| . The generalized T-matrix 
approximation shows an excitation gap around the Fermi energy, which leads to 
similar gaps in the one-particle and two-particle spectral functions. This feature is 
very comforting since it is a manifestation of the energy gap for the excitation of 
particle pairs present in the superfluid. 

This work was partly supported by the KBN under Grant 2P03B02019. 
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Figure 1: The generalized T- matrix equation in the superfluid. 
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Figure 2: The normal part of the self-energy in the T— matrix approximation. 
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Figure 3: The off-diagonal part of the self-energy in the generalized T— matrix 
approximation with the BCS contribution. 
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Figure 4: Diagrams contributing to the generating functional $ in the generalized 
T— matrix approximation. 
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Figure 5: Imaginary part of the diagonal (solid line) and off-diagonal (dashed line) 
part of the generalized T-matrix in the 3 Si channel 
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Figure 6: Two-body spectral function for normal propagators obtained form the gen- 
eralized T-matrix (solid line) and from two uncorrelated BCS propagators (dashed 
line). Projected on the 3 Si channel (upper panel) and 1 Sq channel (lower panel). 
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Figure 7: Two-body spectral function for anomalous propagators obtained form the 
generalized T-matrix (solid line) and from two uncorrelated BCS anomalous prop- 
agators (dashed line). Projected on the 3 S\ channel (upper panel) and 1 So channel 
(lower panel). 
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Figure 8: Two-body spectral function for normal propagators obtained form the gen- 
eralized T-matrix (solid line) and from two uncorrelated BCS propagators (dashed 
line) for nonzero total momentum of the pair (P = 525MeV), projected on the 3 S± 
channel. 
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Figure 9: Two-body spectral function for anomalous propagators obtained form 
the generalized T-matrix (solid line) and from two uncorrelated BCS anomalous 
propagators (dashed line) for nonzero total momentum of the pair (P = 525MeV), 
projected on the 3 Si channel. 
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Figure 10: Real part of the determinant of the inverse T-matrix for the 3 S± (solid 
line) and the 1 S'o channel (dashed line). 
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Figure 11: The position of the quasi-particle pole for normal self-energy for the 
HF propagator (solid line) and for the propagator including the normal self-energy 
from the T-matrix calculation (dashed line). The positions of quasi-particle poles 
of the full propagator including the normal and the off-diagonal self-energies in the 
T-matrix approximation are denoted by the dotted lines. 
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Figure 12: The positions of quasi-particle poles of the full propagator including the 
normal and the off-diagonal self-energies in the T-matrix approximation are denoted 
by the dashed lines. The solid line represents the positions of the poles when taking 
the BCS form only for the off-diagonal self-energy 
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Figure 13: The imaginary part of the retarded self-energy as function of energy for 
p= 0, 140,280,420MeV (solid, dashed, dotted, dash-dotted lines). 
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Figure 14: The single-particle width at the quasi-particle pole of A(p, u) (solid line) 
and at the quasi-particle poles of A s (p, uo) (dashed lines). 
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Figure 15: The real part of the T-matrix contribution to the superfluid gap as 
function of energy for p — 0, 116, 233, 350MeV (solid, dashed, dotted and dash- 
dotted lines). 
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Figure 16: The imaginary part of the T-matrix contribution to the superfluid gap 
as function of energy for p — 0, 116, 233, 350MeV (solid, dashed, dotted and dash- 
dotted lines). 
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Figure 17: The value of the superfluid gap at the quasi-particle pole of u p (dotted 
line) and at the pole E p (dashed line) compared to the BCS gap (solid line). 
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Figure 18: The spectral function A(p,u) including only the diagonal self-energy 
(solid lines) and the full spectral function A s (p, u) (dashed lines) as function of the 
energy. 
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Figure 19: The off-diagonal spectral function B(p,u) for p = 0,280 and 350MeV 
(dotted, solid and dashed lines respectively). 



